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INTRODUCTION 


This  report  is  divided  into  four  sections.  Section  1  is  a  summary  of  progress  and  a  list  of 
project  publications  in  the  past  year.  Section  2  is  a  paper  submitted  to  a  leading 
electrochemistry/corrosion  journal  reporting  on  the  generation  of  hydrogen  in  an  anodicaUy 
polarized  sensitized  ferritic  stainless  steel.  The  paper  demonstrates  that  cavities,  in  this  case 
grooves  along  the  grain  boundary  that  are  created  by  anodic  dissolution,  can  develop  a  critical 
aspect  ratio  that  shifts  the  local  electrode  potential  into  the  range  for  hydrogen  evolution  and 
hydrogen  absorption  into  the  steel.  Section  3  is  a  paper  presented  at  an  international  conference  on 
hydrogen  in  materials  which  treats  the  effects  of  charge  and  mass  transport  limitations  on  the 
distributions  of  cathodic  and  anodic  reactions  on  the  flanks  of  environmentally  induced  cracks. 
Section  4  is  a  reprint  of  a  paper  published  in  a  surface  science  journal  describing  the  atomistics  of 
hydrogen  adsorption  from  the  gas  phase  on  clean  silicon  (100)  surfaces. 

SECTION  1 

PROGRESS  SUMMARY 

The  ongoing  objective  of  this  research  project  is  to  identify  and  evaluate  the  key  factors 
which  affect  the  process  of  hydrogen  entry  into  metals  (particularly  steels)  and  their  subsequent 
embrittlement  under  various  service  conditions,  and  to  develop  adequate  remedial  measures.  This 
includes  (i)  identification  of  the  conditions  which  lead  to  the  generation  of  hydrogen  within  surface 
cavities  that  exist,  for  example,  at  the  junctions  of  metals  with  other  materials  or  at  the  surface 
intersections  of  inclusions  in  metals,  or  that  form  during  service,  e.g.,  crack  propagation  under 
open  circuit  or  cathodic/anodic  protection  situations,  (ii)  the  roles  of  metallurgical  microstructure 
and  solution  species  (H2S)  that  promotes  hydrogen  adsorption  in  duplex  stainless  steels  (AL 
2205),  and  (iii)  the  influence  of  sulfate  reducing  bactering  on  the  cathodic  process  and  absorption 
of  hydrogen  in  steel  (in  collaboration  with  D.  D.  Macdonald’s  ONR  project). 

An  evaluation  of  the  importance  of  the  aspect  ratio  of  cracks  and  other  cavities  has  been  one 
of  the  parameters  under  study  during  the  past  year.  Aspect  ratio  is  defined  as  the  ratio  of  the  depth 


to  opening  dimension  of  a  cavity.  The  aspect  ratio  can  be  independent  of  the  actual  depth  of  a 
cavity.  Hence,  a  very  shallow  crack  and  a  very  deep  crack  may  have  the  same  aspect  ratio.  This 
result  is  readily  obtained  based  on  a  modeling  analysis  in  one  of  our  earlier  ONR  sponsored 
research  publications  (Ateya  and  Rckering,  1975).  Consequently,  a  shallow  crack  could,  in 
principle,  have  the  same  IR  drop  as  a  deep  crack  since  the  magnitude  of  the  IR  is  a  function  of  the 
aspect  ratio  rather  than  simply  of  the  crack  depth.  The  consequences  of  this  independence  of 
aspect  ratio  and  cavity  depth  has  led  to  a  working  hypothesis  on  how  stabilization  of  pit  growth 
following  passive  film  breakdown  may  occur  (Xu  et.  al.,  in  the  list  of  Publications  on  the  Project). 
The  concept  is  equally  powerful  for  the  analysis  of  crack  formation  and  stabilization  during  both 
cathodic  and  anodic  polarization  where  either  hydrogen  absorption  or  metal  dissolution  may  be 
responsible  for  the  cracking  event.  A  discussion  of  the  aspect  ratio  in  the  context  of  potential 
distribution  within  localized  corrosion  processes  is  presented  elsewhere  (Pickering,  J.  Mat.  Sci  & 
Eng.). 

Our  focus  on  the  generation  of  hydrogen  gas  within  growing  surface  cavities  (e.g., 
crevices  and  pits),  and  its  subsequent  absorption  and  embrittlement  of  the  metal  involves 
considerable  amounts  of  mathematical  modeling  and  experimental  testing.  Available  results  have 
been  published  in  several  key  papers  during  the  last  few  years  (Kelly,  et  al.,  Nystrom,  eL  al.,  Xu, 
et.  al.,  Ateya,  et.  al.),  two  of  which  constitute  Sections  2  and  3  of  this  report  and  the  others  are  in 
the  Jan.  93  Annual  Report  and  Final  Report  for  10/1/90  to  9/30/93.  Some  of  the  consequences  of 
these  studies  for  metallic  surfaces  under  cathodic  polarization  impact  on  the  need  for  an  improved 
analytical  characterization  method  for  predicting  (i)  the  success  or  failure  of  cathodic  protection  of 
recesses,  (ii)  the  occurrence  of  active  dissolution  at  the  crack  tip  and  adjacent  crack  flanks,  and  (iii) 
the  distribution  of  hydrogen  charging  rate  on  the  crack  flanks. 
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SECTION  2 


HYDROGEN  GENERATION  IN  AN  ANODICALLY  POLARIZED  SENSITIZED 

FERRITIC  STAINLESS  STEEL 

A.  Sehgal,  B.  G.  Ateya,  R.  N.  Iyer  and  H.  W.  Pickering 

Department  of  Materials  Science  and  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Abstract 

This  paper  presents  evidence  for  the  electrolytic  generation  of  hydrogen  within  the  grain 
boundary  grooves  inside  membranes  of  sensitized  430  stainless  steel  while  the  sample's  outer 
surfaces  are  anodically  polarized  with  respect  to  the  equilibrium  potential  of  the  hydrogen  evolution 
reaction  (HER).  The  grooves  are  formed  by  the  dissolution  of  the  chromium  depleted  zone 
adjacent  to  the  precipitated  chromium  carbides.  This  indicates  that  the  local  electrode  potential 
within  the  grain  boundary  grooves  is  considerably  less  noble  than  that  applied  at  the  outer  surface 
and  indeed  lower  than  that  of  the  HER.  An  IR  potential  shift  mechanism  explains  the  results.  It  is 
shown  that  this  IR  phenomenon  creates  the  conditions  for  hydrogen  embrittlement  within  the  steel 
for  the  existing  electrochemical  conditions  at  the  outer  surface  which  prohibit  hydrogen  generation. 


INTRODUCTION 


Recent  work  on  sensitized  430  stainless  steel  (430  SS)  has  proved  the  existence  of  a  second 
corrosion  process  occurring  during  intergranular  corrosion  under  conditions  of  anodic  polarization 
of  the  steel  surface.  The  process  is  characterized  by  greater  subsurface  dissolution  than 

expected,  specifically  anodic  dissolution  of  the  bulk  alloy  occurred  rather  than  only  of  the  Cr- 
depleted  alloy,  and  by  gas  evolution  from  within  the  resulting  grain  boundary  grooves.  These 
results  were  attributed  to  a  localized  corrosion  process  set  up  by  the  IR  drop  within  the  electrolyte 
of  the  grain  boundary  groove.^  In  this  explanation  the  IR  drop  shifts  the  local  electrode  potential 
within  the  grooves  into  the  active  dissolution  region  of  the  bulk  alloy.^*^  The  observed  gas 
evolution  from  within  the  grain  boundary  grooves  is  also  indicative  of  a  local  electrode  potential 
that  is  considerably  less  noble  than  that  applied  at  the  outer  surface  of  the  specimen  if  the  gas  is 
hydrogen.2  The  evolving  gas  was  deduced  to  be  hydrogen  based  on  similar  experiments  on 
pitting^  and  crevice  corrosion  of  iron.'^  In  the  pitting  investigation,^  the  gas  emanating  from  the 
pits  was  collected  and  found  to  be  hydrogen  by  analysis  in  a  mass  spectrometer. 

The  purpose  of  this  work  is  to  establish  that  the  electrolytically  produced  gas  within  the  grain 
boundary  grooves  is  hydrogen  by  detecting  hydrogen  in  the  metal.  This  was  done  by  testing 
sensitized  430  SS  membranes  in  the  Devanathan-Starchurski  cell^  to  detect  the  hydrogen  generated 
within  the  steel  membranes  while  both  of  its  surfaces  were  maintained  at  anodic  potentials  which 
prohibit  the  HER. 


EXPERIMENTAL 

Membranes  of  type  430  SS  (wt  %  Cr  - 1  6.46,  C  -  0.055,  Mn  -  0.48,  Mo  -  0.024,  N  - 
0.017,  Ni  -  0.32,  P  -  0.035,  S  -  0.004,  Si  -  0.36,  Fe  -  balance)  were  cold  rolled  to  0.037  mm 
thickness  and  solution  annealed  at  1200°C  for  1  hour  and  ice  water  quenched.  Sensitization  was 
induced  by  heating  at  550°C  for  two  hours  followed  by  an  ice  water  quench. 

A  Devanathan-Starchurski  cell^  was  used  to  detect  any  electrolytically  produced  hydrogen 
that  may  have  entered  the  metal  from  within  the  grain  boundary  groove.  An  area  of  0.8  cm^  of 
each  surface  of  the  membrane  was  exposed  to  the  electrolyte.  The  "charging"  side  of  the 
membrane  (in  cell  1)  was  anodically  potentiostated  at  a  potential  of  192  mV  (NHE)  in  deaerated  1 
N  H2SO4  (pH  =  0.36).  A  saturated  calomel  electrode  was  used  as  the  reference  electrode  in  cell  1 
along  with  an  agar  salt  bridge  to  prevent  the  ingress  of  chloride  ions  from  the  calomel  electrode  into 
the  electrolyte. 

The  "exit"  surface  of  the  membrane  (in  cell  2)  was  coated  with  Pd  to  prevent  oxidation  of  the 
underlying  membrane.  An  acid  solution  (1  N  H2SO4  +  3  x  lO"^  M  K3Fe(CN)6,  pH  =  0.36)  was 
used  in  cell  2  instead  of  the  normally  used  basic  solution.  The  potassium  ferricyanide  was  added 
to  the  "exit"  solution  to  detect  iron.  If  the  grain  boundary  grooves  produced  on  the  "charging" 
surface  (cell  1)  continued  to  grow,  eventually  penetrating  the  sample  thickness  and  thus,  arriving  at 
the  passively  maintained  "exit"  surface,  then  the  iron  dissolving  in  the  groove  would  be  detected  as 
Prussian  blue  spots  formed  by  the  reaction  between  ferrous  and  ferricyanide  ions.6  The  "exit" 
surface  was  continuously  monitored  for  prussian  blue  spots  using  a  mirror  placed  in  cell  2  and  a 
high  intensity  light  source.  The  presence  of  potassium  ferricyanide  was  not  found  to  affect  the 
polarization  curve  for  the  sensitized  steel.  A  Hg/Hg2S04  electrode  was  used  to  maintain  the 
potential  at  642  mV(NHE)  for  the  solution  in  cell  2.  This  electrode  potential  was  chosen  so  that 
even  if  the  Pd  layer  did  not  fully  protect  the  steel,  both  the  Cr-depleted  and  bulk  alloy  would  be 
passivated  while  ensuring  complete  oxidation  of  any  hydrogen  atoms  arriving  at  the  "exit"  surface. 
Use  of  potentials  ±  1(X)  mV  of  642  mV  (NHE)  gave  identical  results  during  background  current 
decay.  Hence,  it  was  concluded  that  all  the  hydrogen  atoms  arriving  at  the  "exit"  surface  was 
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oxidized  at  this  potential.  Two  Wenking  75  LB  M  potentiostats  were  used  to  maintain  the 
potentials  on  either  side  of  the  membrane.  Purified  argon  gas  was  used  to  deaerate  the  cell 
compartments.  Further  details  about  the  cell,  equipment  and  procedures  are  reported  elsewhere.'^*^ 


RESULTS  AND  DISCUSSION 

Figure  1  shows  the  measured  anodic  polarization  curve  of  the  sensitized  stainless  steel  in  1  N 
H2SO4.  The  figure  shows  a  secondary  anodic  maximum  at  142  mV  (NHE).  This  has  been 
attributed  to  the  superposition  of  two  anodic  currents,  one  from  the  bulk  alloy  and  the  other  (the 
major  current  density)  from  a  chromium-depleted  zone  adjacent  to  the  grain  boundary  The 
second  maxima  extends  from  10  to  300  mV  NHE.^  The  surface  of  the  sensitized  steel  in  cell  1 
was  maintained  at  192  mV  (NHE)  to  form  the  grain  boundary  (g.b.)  grooves  by  dissolution  of  the 
Cr-depleted  zones  which  lie  along  the  grain  boundaries.  This  potential  is  about  200  mV  more 
noble  than  the  equilibrium  potential  of  the  HER  in  the  IN  H2SO4  solution.  Details  of  the 
metallographic  characterization  of  these  grooves  are  reported  elsewhere.^ 

Figure  2  shows  the  hydrogen  permeation  transient  observed  for  a  0.032  cm  thick  sensitized 
430  SS  anodically  polarized  at  192  mV  (NHE).  Three  different  stages  can  be  identified  in  the 
transient.  In  the  first  stage,  that  lasted  about  900  minutes,  no  permeation  current  could  be  detected. 
In  the  second  stage,  that  lasted  ~  250  minutes,  a  very  small  permeation  current  was  detected.  In 
the  final  stage  the  permeation  current  showed  an  exponential  increase.  During  the  first  and  second 
stages  no  prussian  blue  spots  were  observed  on  the  "exit"  surface.  A  blue  spot  became  visible 
during  the  final  stage  of  the  transient  at  1 170  minutes. 

Figure  3  shows  a  cross-sectional  micrograph  of  a  membrane  polarized  at  192  mV  (NHE)  in 
cell  1.  The  experiment  was  stopped  after  a  permeation  current  of  0.05  |iA  was  recorded  in  cell  2. 

It  is  clear  that  the  grain  boundary  grooves  that  form  on  the  surface  in  cell  1  penetrate  only  part  of 
the  thickness  of  the  membrane.  This  micrograph  is  typical  of  all  of  the  samples  except  that  the 
depth  of  penetration  of  the  g.b.  grooves  varied  from  sample  to  sample,  being  deeper  the  longer  the 
time  of  the  test. 

The  permeation  transient  shown  in  Figure  2  is  much  different  from  that  normally  recorded 
under  conditions  of  cathodic  hydrogen  charging  of  the  external  surface  of  a  membrane,^’^  where 
one  obtains  an  s-shaped  curve  with  a  diffusion  limited  plateau.  The  progressively  increasing  rate 


of  hydrogen  permeation  points  to  the  generation  of  hydrogen  over  more  and  more  area  of  the  g.b. 
grooves  as  the  grooves  increase  in  depth  within  the  membrane. 

The  time  lag  observed  in  region  1  is  the  sum  of  two  components  associated  with  the  various 
processes  forming  the  operating  mechanism,  i.e., 

(i)  The  time  required  to  dissolve  the  Cr-depleted  zone  and  to  generate  grain  boundary 
grooves  of  sufficient  depth  and  aspect  ratio^^  to  set  up  a  large  enough  ohmic  potential 
drop  within  the  grooves  so  that  hydrogen  ion  reduction,  hydrogen  absorption  and  the 
HER  could  take  place. 

(ii)  The  time  taken  for  the  hydrogen  absorbed  in  the  metal  lattice  to  diffuse  from  the  location 
where  it  was  absorbed  (within  the  g.b.  groove,  cell  1)  to  the  Pd  coated  "exit"  surface  of 
the  membrane  where  it  is  oxidized  (cell  2). 

The  above  components  of  the  "time  lag"  cannot  be  separated  easily,  even  though  each  of  them  is 
associated  with  a  distinct  process. 


CONCLUSIONS 

1 .  Detection  of  a  hydrogen  permeation  current  (cell  2)  firom  a  sample  whose  outer  surface  (cell 
1)  was  potentiostated  at  an  anodic  potential  (oxidizing  with  respect  to  hydrogen)  proves  that 
the  local  electrode  potential  inside  the  anodically  produced  grain  boundary  grooves  had 
shifted  more  than  200  mV  in  the  less  noble  direction  into  the  potential  regime  of  the  HER. 
Thus,  hydrogen  ion  reduction  and  its  absorption  into  the  metal  produced  the  permeation 
current  while  at  the  outer  surfaces  (cells  1  and  2)  the  HER  was  prohibited.  The  evolution  of 
gas  bubbles  from  within  the  grooves  shows  that  some  of  the  other  adsorbed  hydrogen  atoms 
combined  to  form  hydrogen  gas  inside  the  grooves. 

2.  The  grooves  form  by  the  dissolution  of  the  chromium  depleted  zone  adjacent  to  the 
precipitated  chromium  carbides  in  the  sensitized  stainless  steel.  For  a  depth-to-opening 
dimension  (aspect  ratio)  of  the  g.b.  groove  that  produced  a  sufficiently  large  IR,  the  HER 
and  hydrogen  absorption  occurred  at  that  depth,  and  then  the  absorbed  hydrogen  diffused  to 
the  outer  surfaces  including  the  "exit"  surface  (cell  2)  where  its  arrival  was  detected  by  the 
oxidation  current  in  cell  2.  The  formation  of  the  grooves  contributed  to  the  total  time  lag  for 
arrival  of  hydrogen  at  the  "exit"  surface. 
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FIGURE  CAPTIONS 


Figure  1.  Measured  anodic  polarization  curve  for  sensitized  Type  430  stainless  steel  in  IN 
H2SO4,  and  a  schematic  of  the  HER  curve.. 


Figure  2.  Measured  hydrogen  permeation  current  at  the  Pd  coated  "exit"  surface  for  a  membrane 
with  both  of  its  surfaces  anodically  polarized  with  respect  to  the  HER. 

Ecelll  =  192  mV  and  Ecell2  =  642  mV  (NHE). 

Figure  3.  Cross  section  optical  micrograph  showing  an  example  of  groove  formation  that 
established  a  sufficiently  large  IR  drop  for  the  occurrence  of  the  HER  inside  the 
groove. 


Figure  1 
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DISTRIBUTIONS  OF  ANODIC  AND  CATHODIC  REACTION  SITES 


DURING  ENVIRONMENTALLY  ASSISTED  CRACKING 


B.  G.  Ateya^  and  H.  W.  Pickering 

Department  of  Materials  Science  and  Engineering 
The  Pennsylvania  State  University 
University  Park,  PA  16802 


Abstract 

This  paper  treats  the  effects  of  charge  and  mass  transport  limitations  on  the  distributions  of 
cathodic  and  anodic  reactions  on  the  flanks  of  environmentally  assisted  cracks.  Under  conditions 
of  cathodic  protection  of  the  metal  surface,  it  is  shown  that  the  presence  of  recesses  or  cracks 
leads  to  significant  changes  in  the  potential  and  ionic  concentrations  at  various  distances  into  the 
cavity.  Consequently,  the  rate  of  the  hydrogen  evolution  reaction  decreases  with  increasing 
distance  into  the  cavity.  Furthermore,  for  active  metals  e.g.,  Fe  or  Ni,  metal  dissolution  occurs 
with  progressively  increasing  rates  at  increasing  distance  into  the  crack.  The  variations  of  these 
cathodic  and  anodic  current  distributions  with  distance  within  the  crack  (or  with  the  time  of 
cracking)  depends  on  the  dimension  of  the  opening  of  the  crack  a,  the  diffusivity  of  the  hydrogen 
ion,  Dh+1  and,  hence,  the  conductivity  of  the  electrolyte,  and  the  polarization  at  the  external 
surface  of  the  metal,  11(0).  Analysis  of  the  migration-diffusion  problem  leads  to  the  development 
of  a  characteristic  depth  given  by 


z  =  [- 


FD 


H+ 

loH 


(Pri(O))  1 

exp — 5 — J 


1/2 


cm 


where  c°^.  is  the  acid  concentration,  and  P  and  io  are  the  charge  transfer  coefficient  and  the 

exchange  current  density  of  the  hydrogen  evolution  reaction,  respectively.  Z  has  a  strong  effect 
on  the  behavior  of  the  system. 


^  Permanent  address:  Chemistry  Department,  Faculty  of  Science,  Cairo  University,  Egypt. 
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The  mechanism  of  environmentally  induced  cracking  of  metals  and  alloys  is  attracting  ever- 
increasing  attention.  Several  conference  proceedings  and  review  articles  have  been  published 
(1-3)  which  review  the  literature  and  project  the  present  state  of  understanding  of  this  issue.  The 
phenomenon  necessarily  involves  complex  interactions  between  the  prevailing  environment  and 
the  electrode  potential  at  the  alloy  surface  with  its  microstructural  features  under  the  effects  of 
stress.  It  is  now  widely  recognized  that  the  chemical  and  electrochemical  conditions  within 
cavities  or  cracks  (and  actually  at  the  tip  of  such  cracks)  in  metals  and  alloys  are  different  from 
those  at  the  external  surface  (4-6).  Over  the  years,  this  recognition  led  to  the  evolution  of  a 
general  belief  that  the  anodic  and  the  cathodic  reactions  are  spatially  separated,  i.e.,  when  the 
anodic  reaction  occurs  within  a  pit,  crevice  or  a  crack,  the  cathodic  reaction  occurs  at  the  outer 
surface  and  vice  versa.  This  assumption  has  been  frequently  used  in  modeling  the  cracking  and 
pitting  processes  (7-9).  It  will  be  shown  below  that  well-separated  anodic  and  cathodic  reactions 
is  only  a  limiting  case  which  is  encountered  in  early  stages  of  crack  growth.  The  general  case  for 
base  metals  involves  distributions  of  both  the  anodic  and  cathodic  (hydrogen  evolution)  reactions, 
with  various  intensities  at  different  locations  on  the  internal  surfaces  of  the  cracks,  along  with  the 
reduction  of  an  oxidant  at  the  outer  surface  for  the  special  case  of  a  freely  corroding  metal.  This 
more  complicated  situation  is  a  result  of  the  combined  effects  of  the  potential  distribution  (IR) 
down  the  crack  and  the  polarization  characteristics  of  the  particular  system.  It  is  also  shown  that 
at  a  later  stage  of  cracking  in  active  (base)  metals,  the  situation  might  well  reach  the  other  limit, 
i.e.,  where  both  the  anodic  and  cathodic  (hydrogen  evolution)  reactions  occur  simultaneously  at 
the  same  place  and  at  the  same  rate  deep  within  the  crack,  under  a  mixed  potential  established  by 
the  particular  metal/hydrogen  evolution  system  existing  inside  the  crack. 

The  objective  of  this  paper  is  to  present  and  explore  the  relations  which  govern  the  distributions  of 
anodic  and  cathodic  reactions  on  the  internal  surfaces  of  cracks  that  propagate  inward  from  the 
sample  surface,  and  the  dependence  of  these  distributions  on  crack  geometry,  the  electrochemical 
parameters  of  the  system,  Ae  properties  of  the  prevailing  environment  and  the  electrode  potential 
at  the  alloy  surface.  These  distributions  are  shown  to  vary  with  the  time  of  propagation  of,  and 
the  distance  within,  the  crack. 

Potential  Prop  within  Cracks:  Causes  and  Consequences 

When  a  cracked  metal  surface  is  held  under  a  certain  electrode  potential,  E(0),  the  flanks  of  the 
crack  are  usually  under  different  "local"  values  of  electrode  potential,  E(x),  where  x  is  the  distance 
into  the  crack.  Figure  1  illustrates  this  situation,  for  a  rectangular  crack  (10).  The  magnitude  of 
the  potential  difference  is  given  by 

AE(x)  =  E(0)-E(x)  =  (p(x)  (1) 

where  (t)(x)  is  the  local  galvani  potential  difference  in  the  electrolyte  which  fills  the  crack,  i.e.,  <j)(x) 
=  0  at  X  =  0.  The  magnitude  and  sign  of  AE(x)  depend  on  the  distance  down  the  crack  x,  the 
crack  opening  a,  the  magnitude  and  sign  of  the  current  density  at  the  outer  surface  (is)  supported 
by  the  potential  at  the  outer  surface,  E(0),  and  the  resistivity  of  the  electrolyte  path  from  x  =  0  to 
X.  This  potential  difference  occurs  in  the  electrolyte  phase.  It  drives  the  ionic  current  through  the 
resistivity  of  the  electrolyte  path  within  the  crack.  This  ionic  current  is  equivalent  to  the  sum  of 
the  faradaic  currents  which  result  from  the  electrochemical  charge  transfer  reactions  occurring  on 
the  flanks  of  the  crack.  The  magnitude  of  this  potential  difference  is  given  by 

(p(x)=  r  I(x)p(x)dx  (2) 

•'O 

where  I(x)  is  the  local  ionic  current  density  flowing  in  the  x  direction  through  the  electrolyte  that 
fills  the  crack,  and  p(x)  is  the  local  resistivity  of  the  electrolyte  path  at  x.  The  relation  between 
E(x)  and  E(0)  is  given  by: 


19 


Figure  1.  Model  of  a  crack  in  which  the  HER  occurs  on  the  outer  surface  (is)  and  on  the 
crack  walls.  The  bulk  electrolyte  is  H'*’  and  Y*  and  the  outer  (x  =  0)  surface 
conditions  are  E  (0)  and  is- 


E(x)  =  E(0) -  r  I(x)p(x)dx  (3) 

*0 

Equation  3  is  written  such  that  the  ionic  current  density  I  has  a  positive  sign  for  anodic  currents 
and  a  negative  sign  for  cathodic  currents.  The  consequences  of  Equation  3  can  be  far-reaching, 
depending  on  the  magnitude  of  the  IR  term.  For  instance  under  the  local  potential  value  of  E(x), 
one  or  the  other  of  the  following  possibilities  occur: 

i.  The  current  density  i(x)  on  the  flanks  of  the  crack,  representing  the  rate  of  the  half  cell 
reaction  occuring  at  x,  becomes  significantly  different  (up  to  an  order  of  magnitude)  from  that 
at  the  external  surface,  is,  under  E(0).  The  local  current  density  i(x)  may  be  greater  or  smaller 
than  is  at  the  outer  surface,  depending  on  the  particulars  of  the  system. 

ii.  For  cracks  in  active  metals,  e.g.,  Fe  or  Ni,  another  electrochemical  reaction  takes  place  if 
favored  by  its  polarization  behavior  and  the  local  value  of  E(x).  In  that  respect,  if  the  crack 
tip  is  under  a  potential  E(x),  it  might  well  encounter  electrochemical  reaction(s)  different  in 
nature  and  rate  than  those  occurring  at  the  outer  surface  of  the  metal. 

We  present  below  some  of  the  governing  equations  and  use  them  to  illustrate  the  distributions  of 
anodic  and  cathodic  reaction  rates  at  different  distances  within  cracks  under  cathodic  polarization 
conditions.  The  effects  of  crack  geometry  (aspect  ratio,  L  cr^  in  Fig.  1),  extent  of  polarization, 
and  the  transpon  properties  of  the  electrolyte  on  these  distributions,  will  be  illustrated  and 
discussed. 

The  Test  System 

The  discussion  is  presented  below  for  a  system  of  iron  in  sulfuric  acid,  for  which  the  computed 
polarization  curve  is  shown  in  Figure  2a.  This  figure  was  constructed  using  the  following  rate 
equations  for  both  the  metal  dissolution  and  the  hydrogen  evolution  (HER)  reactions: 


Figure  2. 


Dlustrarion  of  the  partial  anodic  and  cathodic  currents  for  the  system: 
(a)  Fe/O.IM  Fe2+  in  IN  acid,  (b)  Ni/O.IM  Ni++  in  IN  acid, 

(c)  Cu/O.IM  Cu++  in  IN  acid. 


For 


(4) 
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Fe(s) Fe  (gq)  +2e- 

^Fe  =  ioFe  [exp(anTi(Fe)  /  b)  -  exp(-(l  -  a)nri(Fe)  /  b)] 


(5) 


where  ipe  is  the  metal  dissolution  current  per  unit  of  metal  surface  area,  tKFe)  is  the  polarization 
driving  the  metal  dissolution  reaction,  b  =  RT/F  and  a,  n  and  ipe  are  the  charge  transfer 
coefficient,  the  number  of  electrons  and  the  exchange  current  density  of  this  reaction,  respectively. 
For  the  hydrogen  evolution  reaction  (HER)  i.e., 

2H+  +  2e-  ^  H2  (g)  (6) 

the  rate  is  given  by 

iH  =i^H  [exp(-pTi(H)/b)-exp(-(l-p)tl(H)/b)]  (7) 

where  in,  ioH.  P.  and  ti(H)  refer  to  the  current  density,  exchange  current  density,  charge  transfer 
coefficient  and  the  oveipotential  (polarization),  driving  the  HER.  The  polarization  curves  in 
Figure  2  were  computed  using  typical  values,  i.e.,  ioPe  =  ioH  =  10'^  A  cm'2,  b  =  an  =  0.5  and 
Erev  (Fe/Fe2+)  =  -0.470  V(NHE)  which  corresponds  to  0.1  M  Fe^"*"  ions.  'Die  Erev  (H2/H+)  was 
taJcen  to  be  O.OV  (NHE),  i.e.,  pH  =  0.  With  this  form  of  Equations  5  and  7,  the  cathodic 
polarization  has  a  negative  sign  while  the  anodic  has  a  positive  sign.  Note  that  there  is  a  fairly 
broad  region  of  potential  over  which  reactions  4  and  6  occur  simultaneously  albeit  at  different 
rates.  In  principle,  this  range  extends  from  Erev(Fe/Fe^'^)  to  Erev(H2/H‘'’),  i.e.,  from  -470  mV 
(NHE)  to  0  mV  (NHE).  The  free  corrosion  potential  of  this  system  is  the  potential  at  which  both 
the  anodic  metal  dissolution  and  cathodic  hydrogen  evalution  currents  are  equal  (11). 

The  conresponding  situations  for  Ni  and  Cu  are  shown  in  Figs.  2b  and  c.  For  Ni  the  potential 
range,  over  which  both  nickel  dissolution  and  hydrogen  evolution  occur  simultaneously,  (from 
-280  to  0  mV),  is  much  less  than  in  the  case  of  Fe.  Tlie  case  of  Cu  represents  a  limit  in  that  both 
the  anodic  copper  dissolution  and  cathodic  hydrogen  evalution  reactions  cannot  occur 
simultaneously.  In  fact,  they  are  separated  by  about  350  mV  on  the  potential  coordinate.  These 
figures  will  prove  to  be  most  useful  in  a  discussion  of  the  significance  and  consequences  of 
potential  shifts  within  cracks. 

Transpon  Equations 

We  will  consider  a  model  system  of  a  rectangular  crack  in  an  iron  surface  exposed  to  a  dilute  acid 
solution.  Allowance  is  made  for  one-dimensional  mass  transport  by  diffusion  and  ionic  migration 
in  the  x-direction  into  the  crack.  We  also  use  the  transport  relations  of  dilute  electrolytes.  No 
allowance  will  be  made  here  for  the  additional  effects  of  hydrolysis  (localized  acidification)  or 
complexation  of  metal  ions,  gas  bubbles,  formation  of  salt  film,  etc.  However,  we  do  allow  for 
an  important  factor  which  has  been  frequently  ignored,  i.e.,  the  distribution  of  the  electrochemical 
reactions  on  the  walls  of  the  crack.  We  also  formulate  the  equations  assuming  a  pseudo-steady 
state.  This  assumption  is  justified  in  view  of  the  following  order  of  magnimde  c^culation.  Let  us 
consider  that  the  crack  propagates  1  cm  in  107  s.  Then,  taking  an  average  diffusivity  of  ions 

within  the  electrolyte  of  D  -lO'^  cm^s'^  the  root  mean  square  displacement  during  this  time  is 
Vx^~  V 2x10-5x10’^  ~14  cm.  Thus,  the  rate  of  mass  transfer  (14x10’'^  cm  s*^)  is  much  faster  than 
that  of  crack  propagation  for  values  of  the  crack  propagation  speed  lower  than  1()'7  cm  s*L 

Performing  a  mass  balance,  on  the  various  ionic  species,  on  a  space  element  dx,  as  shown  in 
Figure  1,  one  obtains 


Jh- 


=  -D 


H* 


dcn- 


-I-Ch- 


_f_#1^Ih- 

RTdxJ  F 


dx 


(9) 
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RT  dx 


=  0 


(10) 


dx 


+  2Cpga-> 


F  d(p 
RT^ 


(11) 


where  the  j's  refer  to  the  ionic  fluxes  through  the  electrolyte  in  the  crack  as  a  function  of  distance 
X,  c  and  D  refer  to  the  concentration  and  diffusivity,  respectively,  of  the  various  ionic  sp^ies  and 
<p  is  the  potential  drop  into  the  crack  (Equation  1).  With  this  formulation  of  the  problem,  it 
follows  that  9  >  0  for  anodic  polarizations  and  (p  <  0  for  cathodic  polarizations.  All  the  variables 
in  Equations  9  to  1 1,  (i.e.,  the  c’s,  I's  and  (p)  are  functions  of  distance,  x.  The  concentration 
terms  are  also  subject  to  the  electroneutrality  condition,  i.e., 

Ch*  "t"  ~  ^y~  (12) 

These  fluxes  are  supported  by  the  electron  transfer  reactions  taking  place  on  the  flanks  of  the 
crack.  The  rates  of  these  reactions  are  related  to  the  local  value  of  electrode  potential  and  the 
electrode  kinetics  of  the  particular  reaction  by  the  Butler- Volmer  equation.  Considering  the  space 
element,  dx,  in  Figure  1,  one  obtains 

dlH+=  (^ )  iH  (x)dx 

=  (^)ioH  [exp(-pTi(H)/b)-exp(-(l-p)ri(H)/b)]  dx  (13) 

and 

dlFe2+  =(7)ioFe[exp(anTl(Fe)/b)-exp(-(l-a)nTl(Fe)^)]dx  (14) 

where  Ih'*’.  P  and  T|(H)  are,  respectively,  the  ionic  current  density  resulting  from  the  flow  of  H"*" 
ions  wi^in  the  electrolyte  that  fills  the  crack  in  the  x  direction  (Fig.  1),  the  charge  transfer 
coefficient  of  the  HER  and  the  polarization  driving  the  HER.  Similarly,  lFe2+.  oc  and  T|(Fe)  are 
the  corresponding  parameters  for  the  metal  dissolution  reaction.  The  polarizations,  ri(H)  and 
Ti(Fe),  are  given  by 

ri(H)  =  E(0)  -  tp  (x)  +  0.059  pH  (15) 


Tl(Fe)  =  E(0)-9(x)-Erev(Fe/Fe2+)  (16) 

where  the  Erev  (Fe/Fe2+)  is  the  reversible  potential  of  iron  in  the  particular  medium. 

The  boundary  conditions  of  the  system  are  defined  in  terms  of  ionic  concentrations  and  current  or 
potential  at  the  external  surface.  Thus,  at  x  =  0, 

CFe2+ =  0,  Cy.  =  CH+ =  co  (17) 


and  at  E(0), 

iH  and  iFe  =  f  (polarization  curves)  (18) 

We  now  have  six  equations  (9  to  14)  and  six  variables  (Cj^,  Cy-,  Cpg2-h,  1^+,  lFe2+  and  (j)). 

The  literature  has  no  reference  to  a  complete  solution  to  this  system  of  equations.  Solutions  for 
various  limiting  cases  have  been  report^.  Several  such  solutions  assumed  passive  or  inert  flanks 
of  the  cracks,  and  hence  they  ignor^  the  distribution  of  the  electrochemical  reactions  on  these 
flanks  (7-9).  We  present  below  some  of  the  solutions  which  allow  for  these  distributions,  albeit 
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for  some  limiting  cases.  In  order  to  bring  out  the  physical  significance  of,  and  the  inter¬ 
relationships  among,  the  various  parameters  and  variables  of  the  system,  we  will  discuss  the  cases 
where  analytical  solutions  to  the  equations  were  obtained. 

Cathodic  Polarization 

Under  conditions  of  cathodic  polarization  of  the  external  metal  surface  at  a  potential  E(0)  which 
supports  a  cathodic  current  (is)  at  the  outer  surface,  and  assuming  that  only  the  HER  is  taking 
place  within  the  crack,  the  system  is  governed  only  by  Eqs.  9  and  10.  An  analytical  solution  to 
this  system  is  available  in  the  literature  (10).  The  potenti^,  ionic  concentration  and  current 
distributions  are  given  by  (10): 


cosh[(L-x)/Z] 

cosh[L/Z] 


c(x)  =  c° 


cosh[(L-x)/Z] 

cosh[L/Z] 


i(x)  =  is 


cosh[(L-x)/Z] 
cosh  [L/Z] 


where  L  is  the  crack  depth  and  Z  (in  cm)  is  a  characteristic  depth  given  by 


(19) 

(20) 
(21) 


Z  =  [FD^c°^a/i,f^  (22) 

Note  that  the  characteristic  depth  Z  includes  the  system  parameters  relating  to:  (i)  the  electrode 
kinetics  and  potential  driving  the  HER  (via  is  and  Equation  7),  (ii)  the  width  of  the  crack  opening 
a,  which  determines  the  aspect  ratio  L  0*1,  (iii)  the  concentration  and  the  diffusivity  of  hydrogen 

o 

ions  in,  and  conductivity  of,  the  electrolyte  (via  the  product  Dh+  ch+F).  It  will  be  shown  below 

that  the  characteristic  depth  Z  determines  the  behavior  of  the  system.  Note  also  that  the  current  is  - 
is  related  to  the  potential  E(0)  at  the  external  surface  of  the  metal  through  Equation  7.  Upon 
substituting  for  is  from  Equation  7  (assuming  the  polarization  is  sufficiently  large  to  ignore  the 
reverse  reaction),  one  obtains 


FD 


z  =  [- 


H+ 


loH 


cm 


(23) 


Thus,  a  decrease  in  the  magnitude  of  the  crack  opening  dimension  a  and/or  an  increase  in  the 
magnitude  of  cathodic  polarization  ti(0)  which  drives  the  reaction  at  the  external  surface  (and 
which  has  a  negative  sign)  lead  to  a  decrease  of  the  value  of  Z.  A  similar  characteristic  depth  Z 
has  fi^uently  been  us^  in  studying  the  current  distributions  during  the  electroplating  wiAin 
cracks  or  scratches  (12). 

For  the  sake  of  a  quantitative  calculation,  let  us  consider  the  case  of  an  iron  surface  cathodically 
protected  under  a  potential  of  -  0.5V  (NHE)  in  an  electrolyte  solution  of  pH  =  3  at  25®C.  This 
value  is  comparable  to  the  -0.8V  (Cu/CuS04)  which  is  commonly  used  in  the  practice  of  the 
cathodic  protection  of  steel.  Assuming  the  transfer  coefficient  is  P  =  0.5  and  the  exchange  current 
density,  i©,  of  the  HER  is  10*8  A  cm*A  it  follows  that  the  value  of  is  *  5.5  x  lO'^  A  cm*^.  Using 
values  of  Dh"^  =  lO'^  cm^  s*^  c°  =  lO'^  mole  cm'3,  F  =  96487  C  eq’^  and  a  =  10'2  cm,  Z  =  0.042. 
We  have  calculated  the  current  and  potential  distributions  under  various  conditions  for  several 
combinations  of  L  (from  10'^  to  1  cm)  and  Z  (from  0.004  to  0.04  cm).  The  results  are  presented 
below. 
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Figure  3  illustrates  the  potential  and  (cathodic)  current  distributions  within  cracks  of  different 
depths  ranging  from  L=10  |im  to  1  cm,  for  a  value  of  characteristic  depth,  Z  =  0.04  cm.  The 
distance  coordinate  uses  a  logarithmic  scale  to  accommodate  the  broad  range  of  values  of  crack 
length.  Furthermore,  it  gives  a  measure  of  the  cracking  time.  Assuming  a  constant  value  of  the 
crack  propagation  speed,  e.g.,  10'^  cm  s'^  one  can  easily  calculate  the  cracking  time.  To  that 
extent,  the  calculations  give  the  time  variations  of  the  spatial  potential  and  current  distributions  on 
the  flanks  of  the  crack.  The  pH  within  the  crack  also  changes,  becoming  more  alkaline  with 
distance  into  the  crack  as  given  by  Eq.  20.  If  the  aspect  ratio  is  favorable  for  formation  of  Fe2+ 
ions,  the  pH  will  rise  to  a  limiting  value  which  is  the  equilibrium  pH  of  the  hydrolysis  reaction  as 
discussed  elsewhere  (13). 

Figure  3  reveals  that  as  the  distance  within  the  crack  (or  the  cracking  time)  increases,  the  local 
catiiodic  (HER)  current  density  decreases  while  the  magnitude  of  the  galvani  potential  difference 
in  the  electrolyte,  <))(x).  increases.  Since  4)(x)  is  negative,  Eq.  (1)  quires  that  the  local  electrode 
potential  E(x),  sWfts  in  the  (less  negative)  more  noble  direction  with  increasing  distance  into  the 
crack.  Two  potential  scales  are  us^  on  ^e  right  coordinate  of  Figure  3;  one  for  the  IR  drop  in  the 
crack  electrolyte,  <()(x),  and  the  other  for  the  Ic^  electrode  potential,  E(x).  As  the  distance  within 
the  crack  increases  (for  constant  Z),  the  magnitude  of  ())(x)  increases  (while  (j^x)  is  negative  in 
value),  i.e.,  we  have  a  greater  shift  of  the  local  potentii  E(x)  in  the  more  noble  Oess  negative) 
direction.  No  iron  dissolves  as  long  as  E(x)  is  negative  of  Erev  (Fe/Fe^+X  which  is  -0.470V 
(NHE)  for  the  present  case  while  hydrogen  is  evolved  profusely  under  this  potential.  As  E(x) 
becomes  equal  to  or  positive  of  Erev  (Fe/Fe^'*’),  iron  dissolution  occurs  along  with  hydrogen 
evolution,  at  rates  dependent  on  the  position  of  E(x)  on  the  potential  coordinate  in  Figure  2.  As 
E(x)  shifts  further  towards  less  negative  values,  iH(x)  decreases  while  iFe(x)  increases  gradually. 
At  tiie  limit  of  in  (x)  =  ipe  00  the  value  of  E(x)  equals  the  free  corrosion  potential  of  the  system 
existing  inside  the  crack,  as  shown  in  Figure  2.  One  can  easily  use  Equations  5  and  7  to  calculate 
in  (x)  and  ipefx),  respectively,  at  various  distances  into  the  crack. 

Figure  3  shows  that  under  these  conditions  iron  dissolution  occurs  at  distances  into  the  crack  of  x 
>  0.04  cm  i.e.,  400  |j.m.  Thus,  for  cracks  of  L=10  or  100  pm  depth  under  the  present  conditions 
(i.e.,  a  =  10*2  cm,  E(0)  =  -0.5V  (NHE),  the  entire  depth  of  the  crack  will  be  cathodically  protected 
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Figure  3.  Potential  and  HER  current  distributions  within  cracks  of  various  depths  for  a 
value  of  the  characteristic  depth  Z  =  0.04  cm  during  cathodic  polarization  of  the 
iron  surface  (see  text). 


Local  polenfiol, E( X ),mV (NHE) 
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(see  below  for  a  «  lO'^cm  or  E(0)  <  -0.5V(NHE)).  On  the  other  hand,  for  a  crack  of  0.1  cm 
depth  (1(XX3  |im),  only  the  top  60%  of  its  flanks  will  be  under  cathodic  protection.  The  bottom 
40%  will  be  undergoing  simultaneous  anodic  dissolution  and  hydrogen  evolution  at  rates  which 
are  dependent  on  distance,  as  discussed  above.  Alternatively,  for  a  1  cm  deep  crack,  only  the  top 
4%  of  its  depth  is  under  cathodic  protection  while  the  rest  is  undergoing  ano^c  dissolution. 

Using  the  E(x)  value  calculated  from  Eq.  1  and  a  knowledge  of  the  active/passive  behavior  as  a 
function  of  pH,  it  can  be  determined  if  the  anodic  metal  dissolution  undergoes  active  or  passive 
dissolution  as  a  function  of  distance  into  the  crack  (or  cracking  time),  as  reported  elsewhere  for 
crevice  corrosion  (14). 

Figures  4  and  5  illustrate  the  effects  of  one  order  of  magnitude  decrease  in  Z,  i.e.,  for  Z  =  0.004 
cm.  This  decrease  in  Z  was  brought  about  in  two  different  ways,  i.e.,  by  decreasing  a  from  lO'^ 
to  10^  cm  in  Fig.  4,  keeping  E(0)  =  -0.500V  (NHE)  and  all  other  parameters  constant.  Figure  4; 
and  by  changing  E(0)  from  -0.50  to  -0.72  (NHE)  while  keeping  a  =  lO'^  cm  and  the  other 
parameters  constant.  Figure  5.  For  the  case  of  a  =  10^  cm  in  Fig.  4,  the  distance  beyond  which 
iron  dissolution  occurs  becomes  much  shoner,  being  in  this  case  only  about  0.(X)4  cm,  i.e.  40 
|im.  Under  this  condition  iron  dissolution  occurs  at  distances  greater  than  40  ^.m.  Consequently, 
for  any  crack  of  a  depth  less  than  40  |im  (for  the  conditions  of  Fig.  4),  the  entire  crack  depth  will 
be  under  cathodic  protection.  Similarly,  for  the  conditions  in  Fig.  5,  cracks  will  be  fully  protected 
if  less  than  400  |im  in  depth.  Conversely,  for  a  1  cm  deep  crack,  only  the  top  4%  of  its  depth  will 
be  under  cathodic  protection  while  its  interior  (96%  of  its  depth)  will  be  involved  in  anodic 
dissolution. 

For  very  narrow  cracks  (<3=100  A  and  E(0)  =  -0.5  (NHE))  only  the  top  1.6  |im  of  cracks  of  any 
depth  will  be  under  cathodic  protection,  the  rest  being  involved  in  anodic  dissolution. 

Conclusions 


The  present  analysis  of  the  diffusion-migration  problem  within  environmentally  assisted  cracks 
leads  to  the  following  conclusions: 


1. 


The  current  and  potential  are  nonuniformly  distributed  within  cracks.  A  characteristic  depth 
Z  was  derived  which  combines  the  important  parameters  of  the  system,  and  affects  these 
distributions. 


Z=[ 


k)H 


(pri(0))-,l^ 

exp  — r —  J  cm 


2.  For  conditions  of  cathodic  polarization  at  the  external  metal  surface,  the  potential  within  the 
crack  shifts  in  the  noble  direction  towards  a  limiting  value,  as  a  result  of  the  potential  shift 
(IR)  mechanism.  The  nature  of  this  limiting  value  of  potential  is  different  for  active  than  for 
noble  metals. 


3 .  This  potential  shift  leads  to  significant  decreases  in  the  rate  of  hydrogen  evolution  on  the 
flanks  of  the  crack. 


4 .  For  active  metals,  the  potential  shift  may  become  sufficiently  large  to  promote  anodic 
dissolution  within  the  crack  even  though  the  external  metal  surface  is  under  cathodic 
protection. 

5 .  The  most  general  case  involves  nonuniform  distributions  of  both  the  cathodic  and  anodic 
reactions  on  the  flanks  of  cracks. 


6.  An  initial  crack  in  an  active  metal  may  remain  under  cathodic  protection  until  its  aspect  ratio, 
L  or  depth  L  (for  a  constant  opening  dimension  a)  reaches  a  certain  critical  value  when 
the  potential  at  L  shifts  to  the  region  of  metal  dissolution.  This  critical  depth  depends  on  the 
value  of  a,  the  properties  of  the  environment  and  on  the  polarization  at  the  metal  surface 
through  the  characteristic  depth  Z. 


Dimensionless  curreni  iu+  ( x  )  il'  Dimensionless  current  i„+  (  x )  i^ 
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Dislonce  {depthi  within  crock, x  =  vt, cm 


Effect  of  the  decrease  in  crack  opening  dimesion  a  on  the  potential  and  current 
distributions  for  the  conditions  in  Fig.  3. 


Disfonce  (depth)  within  crock, x-vt, cm 


Figure  5. 


Effect  of  the  decrease  in  cathodic  polarization  at  the  outer  surface,  E(0),  on  the 
potential  and  current  distributions  for  the  conditions  in  Fig.  3. 


Local  potential, E(  X  ), mV  (NHE)  Local  potential, E(  x), mV  (NHE) 
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A  field  ion-scanning  tunneling  microscopy  study  on  hydrogen  chemisorption  on  the  Si(100X2  x  1)  surface  is  presented.  At  low 
coverages,  hydrogen  atoms  reside  singly  on  top  of  the  dimerised  Si  atoms,  and  are  imaged  brightly.  The  hydrogen  chemisorption 
induces  buckling  of  dimers,  indicating  the  strong  bonding  between  Si  and  H  atoms.  With  increasing  coverage,  both  the  (2x1) 
monohydride  and  (1  x  1)  dihydride  phases  were  formed.  The  former  is  imaged  dark  compared  with  the  unreacted  Si  dimers,  due  to 
the  reduction  of  the  density  of  electronic  states  near  the  Fermi  level.  Surface  etching  was  also  observed.  It  was  found  that  the 
corrosion  of  the  surface  is  modest  in  the  monohydride  phase,  while  during  the  formation  of  the  dihydride  phase,  the  corrosion 
becomes  significant.  The  behavior  of  hydrogen  desorption  from  the  dihydride  and  monohydride  phases  was  investigated  as  a 
function  of  annealing  temperature.  Our  STM  results  support  the  mechanism  that  the  desorbing  H,  molecules  are  formed  by 
combination  of  two  hydrogen  atoms  forming  the  dihydride  phase.  Upon  annealing  at  elevated  temperatures,  the  Si  overlayer  stripes 
are  formed  by  desorption  of  hydrogen  from  the  etching  products  and  the  rearrangement  of  the  Si  atoms. 


1.  Introduction 

The  chemisorption  of  atomic  hydrogen  on  vari¬ 
ous  silicon  surfaces  has  been  studied  for  several 
decades  [1-4]  because  of  its  importance  in  both 
fundamental  research  and  technological  applica¬ 
tions.  For  instance,  it  is  well  known  that  hydrogen 
passivates  silicon  surfaces  by  saturating  the  sili¬ 
con  dangling  bonds,  other  electronically  active 
defects  and  stabilizing  grain  boundaries,  etc.  In 
addition  to  that,  the  adsorption  of  hydrogen  on  Si 
surfaces  also  provides  a  model  example  for  inves¬ 
tigation  of  the  surface  reconstructions  and  bond¬ 
ing  properties.  Though  atomic  hydrogen  adsorp¬ 
tion  on  the  Si(  100X2  x  1)  surface  has  been  exten¬ 
sively  studied  by  various  kinds  of  methods,  many 
controversies  still  exist.  Previous  studies  have 
confirmed  the  existence  of  three  reconstructions 
induced  by  hydrogen  adsorption  at  different  con¬ 
ditions,  i.e.  (2x1)  [2,5]  (the  monohydride  phase), 
(1  X  1)  [2]  (the  dihydride  phase),  and  the  (3  x  1) 
[3]  phase.  By  modest  adsorption  of  atomic  hydro¬ 
gen  at  room  temperature  or  saturate  adsorption 


at  elevated  temperature  (about  250°C),  the  (2x1) 
phase  can  be  formed,  which  is  characterized  by 
binding  one  hydrogen  atom  to  the  dangling  bond 
on  each  dimerized  silicon  atom  without  breaking 
the  dimer  bond.  By  saturating  the  surface  with 
hydrogen  at  room  temperature,  Sakurai  and 
Hagstrum  [2]  discovered  the  (1x1)  reconstruc¬ 
tion  and  introduced  the  idea  of  the  dihydride 
phase,  which  formed  by  adsorbing  two  hydrogen 
atoms  on  each  Si  atom  accompanied  by  the  inter¬ 
ruption  of  the  dimer  bond.  Yet,  the  geometry  of 
the  hydrogen  atoms  in  this  phase  has  not  been 
identified  unambiguously.  Recently,  Chabal  and 
Raghavachari  [3]  found  a  new  (3x1)  phase  at 
saturation  coverage  while  holding  the  substrate 
temperature  at  about  100°C.  This  reconstruction 
is  formed  by  an  alternating  arrangement  of 
monohydride  and  dihydride  rows  and  appears 
even  more  stable  than  either  the  monohydride 
phase  or  the  dihydride  phase  over  a  certain  range 
of  chemical  potential  [6].  The  electronic  struc¬ 
tures  of  the  (2  X  1)  monohydride  phase  and  (1  x 
1)  dihydride  phase  have  been  investigated  by 
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Sakurai  and  Hagstrum  [2]  and  Ciraci  et  al.  [7] 
using  UPS,  and  by  Maruno  et  al.  [8,9]  using 
EELS,  etc.  It  was  found  that  the  hydrogen  in¬ 
duced  electronic  states  are  located  at  about  — 12 
and  - 10  eV  with  respect  to  the  vacuum  level  for 
the  monohydride  phase  and  dihydride  phase,  re¬ 
spectively  [2].  From  TPD  (temperature  pro¬ 
grammed  desorption)  [10]  experiment,  Gupta  et 
al.  derived  that  desorption  of  Hj  from  both  the 
monohydride  and  dihydride  phases  exhibits  sec¬ 
ond-order  kinetics.  However  a  previous  investiga¬ 
tion  [11]  and  more  recent  LITD  (laser  induced 
thermal  desorption)  experiment  [12]  suggested  a 
first-order  desorption  kinetics  for  desorption  of 
H2  from  both  phases. 

The  corrosion  of  the  surface  upon  hydrogen 
adsorption  has  also  been  noticed  and  investigated 
mainly  by  desorption  experiments.  The  products 
of  the  reaction  have  been  identified  to  be  mainly 
SiHj  and  small  fractions  of  SijH^,  SijHg  and 
Si4H]o,  etc.  [13],  but  the  detailed  process  and 
etching  mechanism  are  still  unclear.  Moreover, 
earlier  conclusions  were  based  on  the  data  aver¬ 
aged  over  the  macroscopic  area.  The  observation 
of  the  surface  structure  in  real  space  had  not 
been  possible  until  the  invention  of  the  scanning 
tunneling  microscope  by  Binnig  and  Rohrer  [14]. 
The  STM  has  proven  to  be  of  great  potential  to 
study  the  surface  geometry  as  well  as  the  surface 
electronic  states  with  atomic  resolution  in  real 
space.  Hamers  et  al.  [15]  imaged  the  monohy¬ 
dride  phase  produced  by  dissociative  adsorption 
of  NHj  on  the  Si(l(K)X2  x  1)  surface.  The  imaged 
feature  was  attributed  to  the  tunneling  current 
through  the  localized  Si-H  bonding  orbitals  with 
a  sample  bias  of  —2  V.  However,  Johansson  et  al. 
[16]  concluded,  based  on  their  polarization-de- 
pendent  ARPES  experiment,  that  the  highest  en¬ 
ergy  position  of  this  band  was  4.7  eV  below  the 
Fermi  level  and  that  no  H-induced  structures 
should  be  visible  at  all  in  the  energy  ranging  from 
0  to  -2  eV,  in  contrary  to  Hamers  et  al.  their 
STM  study.  Recently,  Boland  reported  the  obser¬ 
vation  of  the  newly  discovered  (3x1)  phase  [17] 
and  discussed  the  role  of  pairing  in  the  recombi- 
native  desorption  of  hydrogen  [18].  However,  the 
understanding  of  this  system  is  far  from  complete 
and  many  controversies  remain  unresolved. 


In  this  paper,  we  will  present  our  FI-STM 
investigations  on  the  hydrogen  chemisorption  on 
the  Si(100X2  x  1)  clean  surface.  In  the  following 
sections,  after  a  brief  description  of  our  experi¬ 
mental  arrangement,  we  will  present  the  behavior 
of  adsorbed  hydrogen  atoms  at  the  initial  stage  of 
adsorption,  then  we  will  show  the  surface  struc¬ 
ture  of  the  monohydride  phase  and  dihydride 
phase,  and  the  corrosion  of  the  surface  caused  by 
hydrogen  adsorption.  Finally  we  will  show  the 
surface  change  upon  desorption  at  different  tem¬ 
peratures.  The  mechanism  behind  the  presented 
results  will  also  be  discussed. 


2.  Experiment 

The  details  of  our  Fl-STM  system  have  been 
described  elsewhere  [19].  Here  is  a  brief  descrip¬ 
tion.  Our  Fl-STM  is  a  combination  of  a  high 
performance  scanning  tunneling  microscope 
(STM)  and  a  room  temperature  field  ion  micro¬ 
scope  (FIM)  which  is  used  for  the  precise  charac¬ 
terization  of  the  scanning  tip  and  its  atomic  scale 
fabrication.  The  entire  STM  setup  is  suspended 
with  three  1(X)  cm  long  springs  and  a  12-piece  set 
of  magnetic  eddy  current  dampers  which  is  at¬ 
tached  beneath  the  STM  set.  The  entire  Fl-STM 
chamber  is  mounted  on  the  so-called  active 
damper.  By  doing  so,  the  troublesome  resonance 
frequencies  at  1-3  Hz  are  completely  eliminated, 
making  it  possible  to  operate  the  STM  stably  and 
reliably.  The  use  of  the  (111)  oriented  single 
crystal  tungsten  tip  provides  us  with  needed  sta¬ 
bility  and  long-lifetime  of  several  months  with 
daily  operation.  The  main  UHV  chamber  is 
pumped  by  a  1200  (/i  custom-design  Riber  ion 
pump  and  the  base  pressure  is  below  3  x  10“  ” 
Torr  except  for  the  time  of  annealing  the  sample 
up  to  1220‘’C.  The  pressure  rise  during  the  final 
annealing  does  not  go  beyond  1  x  10“’°  Torr 
with  this  fast  pumping  system. 

The  Si(lOO)  sample  used  was  cut  into  a  4  x  19 
mm^  piece  from  a  Sb  doped  n-type  commercial 
SKKX))  wafer  with  p  about  0.07  fl  •  cm.  Prior  to 
introducing  the  sample  into  the  UHV  FI-STM 
chamber,  only  washing  with  methyl-alcohol  was 
administered  for  sample  cleaning.  Cleaning  the 
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sample  in  the  UHV  preparation  chamber  (which 
is  attached  to  the  main  chamber)  was  carried  out 
by  a  series  of  heating  up  to  1220°C  following  the 
5  hour’s  degassing  at  700°C.  In  this  way,  the  clean 
Si(l(X)X2  X  1)  surface  with  a  minimum  density  of 
defects  was  obtained.  The  temperature  was  meas¬ 
ured  by  an  infrared  pyrometer. 

The  atomic  hydrogen  dosing  was  achieved  by 
introducing  high  purity  Hj  gas  into  the  UHV 
chamber  through  an  adjustable  leak  valve,  pass¬ 
ing  through  the  hot  tungsten  ribbon  which  was 
heated  up  to  1500°C  to  dissociate  molecular  hy¬ 
drogen  into  atomic  hydrogen.  The  Hj  dosage, 
which  was  given  in  Langmuir  units,  was  used  only 
as  a  measure  for  the  atomic  hydrogen  deposition 
rate.  It  was  ensured  that  the  hot  tungsten  ribbon 
did  not  damage  the  vacuum  during  H  exposure. 
The  hydrogen  adsorption  was  performed  by  hold¬ 
ing  the  sample  at  room  temperature.  Some  tem¬ 
perature  rise  in  the  sample  surface  was  expected 
during  the  H  adsorption,  however,  the  tempera¬ 
ture  rise  was  not  high  enough  to  cause  desorption 
of  atomic  hydrogen  from  the  surface.  The  STM 
images  were  taken  at  least  3  h  after  the  annealing 
to  avoid  the  effect  of  thermal  drift.  A  sample  bias 
between  -3  and  -8  3  V  was  used,  the  tunneling 
current  was  kept  constant  at  20  pA  and  the 
scanning  rate  was  450  A/s  during  the  normal 
STM  operation. 


3.  Results  and  discussion 

3.1.  Properties  of  the  initial  stage  of  hydrogen  ad¬ 
sorption 

Fig.  1  is  a  typical  STM  topography  of  atomic 
hydrogen  adsorbed  on  the  SKI 00X2  x  1)  surface 
at  a  very  low  coverage  ( ~  0.03  ML).  The  surface 
is  characterized  by  the  bright  round  spots  sitting 
on  top  of  the  buckled  dimers  and  distributed 
randomly  on  the  surface.  At  both  sides  of  these 
bright  spots,  the  dimers  are  also  buckled,  imaged 
brighter  than  the  unbuckled  dimer  chains.  Based 
on  the  fact  that  the  number  density  of  these 
bright  spots  increases  with  the  hydrogen  dosage, 
these  are  concluded  to  be  the  image  of  hydrogen 
atoms.  This  is  basically  consistent  with  the  result 


Fig.  1.  STM  topography  of  Si(100X2xl)  surface  covered  by 
about  0.03  ML  of  atomic  H,  sample  bias  K,  =  1.6  V. 


reported  by  Boland  [18],  though  he  could  not 
observe  the  buckling  of  the  dimers  adjacent  to 
the  hydrogen  atoms.  A  question  is,  is  it  hydrogen 
atoms  that  prefer  to  stay  on  buckled  dimers  or  is 
the  buckling  of  the  dimers  induced  by  the  ad¬ 
sorbed  hydrogen  atoms?  From  the  following  ar¬ 
gument,  we  conclude  that  it  is  hydrogen  atoms 
that  induce  the  buckling  of  the  dimers  adjacent  to 
them,  though  it  is  still  possible  that  H  atoms 
prefer  to  stay  on  buckled  dimers.  First,  by  com¬ 
paring  the  density  of  the  buckled  dimers  before 
and  after  adsorption,  it  was  found  that  the  den¬ 
sity  of  this  kind  of  buckled  dimers  increased  with 
increasing  H  dosage.  Second,  adatom-like  protru¬ 
sions  without  buckling  of  the  adjacent  dimers 
were  rarely  found.  The  buckling  of  the  dimers  on 
the  Sid 00X2  X  1)  surface  has  been  of  interest  for 
many  years  and  there  are  controversies  between 
the  results  obtained  from  STM  and  other  tech¬ 
niques.  Results  from  the  STM  [20]  have  sug¬ 
gested  that  the  surface  seems  mostly  unbuckled. 
However,  it  has  been  revealed  that  the  energy 
difference  between  the  buckled  dimer  and  the 
unbuckled  dimer  is  small  enough  so  that  the  STM 
imaged  symmetric  dimers  may  in  fact  be  the  time 
average  of  rapidly  oscillating  asymmetric  (buck¬ 
led)  dimers  with  charges  transferring  back  and 
forth  between  two  dimerized  atoms  [21,22],  while 
dimers  can  be  stably  buckled  in  the  vicinity  of 
steps  and  various  kinds  of  point  defects  [20].  The 
fact  that  the  H  atoms  are  almost  adsorbed  singly 
instead  of  in  pairs  at  low  coverage  indicates  the 
significant  strong  bonding  between  the  H  atom 
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and  the  dangling  bond  of  one  of  the  dimerized  Si 
atoms.  This  also  indicates  that  most  of  the  arrived 
H  adatoms  adsorbed  on  the  dangling  bonds  with¬ 
out  moving  around.  So  the  periodicity  and  the 
symmetry  around  the  adsorption  site  no  longer 
exist  due  to  this  strong  bonding  as  well  as  the 
saturation  of  the  dangling  bond  by  adsorbed  hy¬ 
drogen  atom.  This  causes  a  charge  redistribution, 
which  is  likely  to  induce  the  buckling  of  the 
dimers  nearby.  Similar  results  were  also  obtained 
for  alkali-metal  adsorption  on  the  Si(100X2  x  1) 
surface  [23]  by  STM.  From  the  number  of  buck¬ 
ling  dimers  induced  by  H  adsorption,  we  suggest 
that  the  interaction  between  the  H  atom  and  the 
Si  substrate  is  of  short  range,  i.e.,  limited  to  one 
or  two  dimers  near  the  adsorbed  H  atom.  It  is 
interesting  and  important  to  study  the  interaction 
between  the  singly  adsorbed  hydrogen  atoms  with 
surface  silicon  atoms  as  well  as  the  interactions 
among  the  adsorbed  hydrogen  atoms  through  the 


investigation  of  the  initial  stage  of  hydrogen  ad¬ 
sorption,  which  may  reveal  the  mechanism  of  the 
desorption  process  [18].  However,  up  to  now, 
little  work  has  been  done  on  it,  which  limits  our 
further  discussion.  We  hope  our  results  will  stim¬ 
ulate  further  investigations  of  the  detailed  prop¬ 
erties  of  the  Si-H  bond,  such  as  the  charge 
transfer,  bond  length,  electronic  states,  etc. 

3.2.  Formation  of  the  monohydride  phase  and  di¬ 
hydride  phase 

Figs.  2a-2d  are  the  STM  images  showing  the 
changes  of  the  surface  structure  with  increasing 
Hj  coverage.  Fig.  2a  is  the  image  obtained  after 
24  L  adsorption,  the  coverage  is  estimated  to  be 
approximately  0.3  ML.  In  this  case,  the  surface  is 
characterized  by  the  interruption  of  the  dimer 
chain  by  dark  areas  and  the  buckling  of  almost  all 
the  dimers  on  the  surface  which  imaged  bright. 


Fig.  2.  Si(lOO)  surface  structures  after  hydrogen  adsorption  at  different  dosages,  (a)  After  adsorption  of  24  L,  K,  =  - 2.0  V,  (b)  after 
adsorption  of  48  L,  K,  -  -3.0  V,  (c)  after  adsorption  of  108  L,  -  1.6  and  (d^)  after  adsorption  of  2000  L,  -2.0  V.  The 
area  is  about  100  x  100  A"  for  (a)-(c)  and  about  43  x  43  A^  for  (d),  respectively. 
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Fig.  3.  The  illustration  of  two  different  kinds  of  features, 
marked  with  “A”  and  “B”,  respectively.  The  image  was 
obtained  upon  48  L  exposure.  =  -2.0  V. 


With  increasing  coverage,  the  area  of  dark  con¬ 
trast  increases,  as  shown  in  fig.  2b,  which  was 
obtained  upon  48  L  exposure.  In  this  case,  the 
dimer  chains  in  the  dark  area  are  well  resolved, 
showing  clearly  that  the  dimers  in  the  monohy¬ 
dride  phase  are  imaged  dark.  This  can  be  ex¬ 
plained  by  the  reduction  of  the  electronic  density 
of  states  near  the  Fermi  level  [25,26].  From  the 
UPS  results  [7,24]  and  theoretical  calculation  [27], 
it  is  known  that  on  the  clean  dimerized  Si(lOO)- 
(2x1)  surface,  there  exists  a  localized  surface 
state  near  the  Fermi  level  due  to  the  dangling 
bond  of  the  surface  silicon  atoms.  This  surface 
state  is  removed  upon  hydrogen  adsorption  by 
saturation  of  the  dangling  bonds,  resulting  in  the 
reduction  of  the  tunneling  current. 

When  the  surface  is  largely  covered  by  the 
monohydride  phase,  the  prominent  feature  is  the 
bright  dots  on  the  well  ordered  (2x1)  monohy¬ 
dride  surface  (see  also  fig.  2b).  These  bright  dots 
can  be  grouped  into  two  types,  namely  type  A 
and  B.  Fig.  3  is  a  close-up  illustrating  typical 
images  of  dots  “A”  and  “B”,  where  “A”  dots  are 
located  on  the  side  of  the  dimer  chains,  while 
“B”  dots  are  located  on  the  dimer  chains  and 
imaged  larger  than  “A”  dots.  According  to  the 
observation  that  (a)  the  dimers  adjacent  to  “A” 
dots  are  rarely  buckled,  which  is  in  contrast  to 
the  results  in  the  initial  stage  of  hydrogen  adsorp¬ 
tion,  and  that  (b)  the  density  of  “A”  dots  de¬ 
creased  with  the  increase  of  the  hydrogen  cover¬ 
age,  we  suggest  that  “A”  dots  represent  the  un¬ 


paired  hydrogen  atoms.  Thus  the  unbuckling  of 
the  dimers  adjacent  to  the  “A”  dots  indicate  the 
considerable  passivation  of  the  surface  reactivity. 
On  the  other  hand,  the  density  of  type  “B”  dots 
remains  almost  unchanged.  Further  increasing  the 
dosage  results  in  the  formation  of  small  amounts 
of  clusters  imaged  slightly  larger  than  “B”  dots. 
Though  the  identification  of  the  products  by  STM 
alone  is  not  possible,  we  suggest  that  these  type 
“B”  dots  are  the  products  of  reacted  Si  atoms 
with  hydrogen,  according  to  the  thermal  desorp¬ 
tion  results  by  Gates  et  al.  [13],  mainly  SiHj  with 
small  amount  of  other  species,  such  as  Si2H6, 
SijHg  and  Si4H,o.  Fig.  2c,  was  obtained  upon 
108  L  of  hydrogen  exposure,  where  type  “A”  dots 
were  almost  eliminated,  while  type  “B”  dots  re¬ 
mained.  It  is  also  found  from  our  STM  investiga¬ 
tion  that  these  products  spread  randomly  and 
almost  separately  on  the  surface,  with  little  ten¬ 
dency  to  combine  into  larger  clusters,  indicating 
the  absence  of  etching  centers  or  nucleation  cen¬ 
ters  of  the  reacted  Si  atoms.  During  the  forma¬ 
tion  of  the  monohydride  phase,  the  degree  of 
surface  etching  is  still  limited  within  a  low  level, 
which  let  the  monohydride  phase  form  with  higher 
order.  Also  in  the  case  as  shown  in  fig.  2c,  it  was 
found  that  the  formation  of  a  new  phase,  the 
dihydride  phase,  started,  indicated  by  the  split¬ 
ting  of  the  monohydride  dimers. 

With  increasing  dosage,  the  surface  gradually 
shifts  from  the  (2  x  1)  monohydride  phase  to  the 
dihydride  phase.  Fig.  2d  shows  the  surface  upon 
2000  L  of  hydrogen  adsorption.  The  surface  is  of 
the  (1  X  1)  reconstruction  with  a  fraction  of  the 
(2  X  1)  phase.  We  noticed  that  with  the  increase 
of  the  proportion  of  the  (1  x  1)  phase,  the  etching 
of  the  surface  becomes  more  and  more  extensive, 
shown  by  the  high  density  of  missing  atom  defects 
and  “B”  dots.  The  difference  between  the  sizes 
of  the  dots  is  slight,  indicating  that  the  bonding  to 
the  substrate  is  still  strong  enough  to  resist  the 
migration  of  the  reacted  species. 

3.3.  Desorption  of  hydrogen 

The  understanding  of  the  desorption  proper¬ 
ties  of  hydrogen  is  almost  of  the  same  importance 
as  the  adsorption  process,  especially  in  control- 
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ling  of  the  epitaxial  growth  of  the  Si  film  on  the 
Si(lOO)  surface  by  CVD,  as  well  as  analyzing  the 
surface  etching  process.  For  many  years,  several 
groups  have  studied  the  properties  of  the  desorp¬ 
tion  process  [10,12,13].  The  desorption  tempera¬ 
tures  for  hydrogen  in  the  SiH  species  (^j-state) 
and  SiH  2  species  Oj-slale)  have  been  deter¬ 
mined  by  several  reports  with  similar  conclusions. 
Gupta  et  al.  [10]  reported  that  for  hydrogen  ad¬ 
sorbed  on  porous  silicon,  which  is  similar  to  that 
on  the  Sid  00)  surface  according  to  the  similari¬ 
ties  between  the  infrared  spectra  of  hydrogen  in 
both  cases,  the  desorption  of  hydrogen  in  the 
/3, -state  occurs  at  temperatures  between  720  and 
800  K  (447  and  527°C),  while  the  desorption  of 
hydrogen  in  the  /Sj-state  occurs  at  640-700  K 
(367-427°C).  From  TPD  spectra.  Gates  et  al.  [13] 
found  that  the  peak  temperature  for  desorp¬ 
tion  from  the  /Sj-state  is  540°C  and  the  desorp¬ 
tion  starts  at  about  450‘’C,  while  the  peak  temper¬ 
ature  for  desorption  of  the  /32-state  is  425°C  and 
the  desorption  starts  at  about  370°C.  So  these 
results  suggest  the  desorption  of  hydrogen  starts 
about  370°C  for  the  ^2-state  and  450'C  for  the 
^,-state.  Gates  et  al.  [13]  also  found  that  the 
desorption  of  the  SiH  3  species  occurs  at  about 
375°C  in  the  form  of  SiH  4,  which  is  50'C  below 
the  desorption  peak  of  the  desorp¬ 

tion  of  SiH  4  also  results  in  the  removal  of  the 
silicon  atoms  from  the  top  of  the  surface.  The 
TPD  yield  of  SiH  4  saturates  at  the  maximum 
value  when  the  coverage  increased  beyond  0.5 


ML.  Recently,  Boland  [18]  reported  STM  data  on 
the  behavior  of  the  monohydride  phase  upon  H 
desorption.  However,  controversies  about  the 
desorption  kinetics  and  transition  process  from 
the  dihydride  phase  to  the  monohydride  phase 
remain.  Information  on  the  real  space  observa¬ 
tion  of  the  surface  structure  change  upon  H2 
desorption  is  still  not  sufficient,  especially  in  the 
process  of  the  desorption  from  the  dihydride 
phase. 

To  investigate  the  desorption  process  of  hydro¬ 
gen  from  both  the  SiH 2  state  {^2^  state 

(Bi),  the  surface  was  exposed  to  1000  L  of  hydro¬ 
gen  prior  to  annealing  at  modest  temperatures. 
As  shown  in  fig.  4a,  the  surface  is  of  mixture  of 
the  (1  X  1)  and  (2  X  1)  phases.  According  to  the 
temperature  dependence  of  the  H2  desorption  in 
each  species,  we  performed  sequential  anneal¬ 
ings,  i.e.  heating  the  sample  for  1  min  at  tempera¬ 
tures:  320“C,  400‘’C,  460°C  and  500“C. 

Annealing  the  sample  at  320°C  resulted  in  no 
obvious  changes  in  the  surface  structure,  which  is 
consistent  with  the  TPD  result  at  the  same  tem¬ 
perature,  that  is,  almost  all  the  adsorbed  hydro¬ 
gen  atoms  and  reacted  species  remain  on  the 
surface.  Though  some  of  the  hydrogen  atoms  may 
obtain  sufficient  energy  to  desorb  from  the  sur¬ 
face  due  to  the  Boltzmann  distribution  of  thermal 
energy,  the  amount  is  small  in  such  a  short  time 
duration  as  was  shown  by  an  isothermal  desorp¬ 
tion  experiment  [10].  After  further  annealing  of 
the  sample  at  400'’C  for  1  min,  as  shown  in  fig. 


Fig  4  Changes  in  surface  structure  upon  annealing  at  different  temperature,  (a)  Mixture  of  monohydrogen  and  dihydrogen  after 
adsorption  of  1000  L  hydrogen  without  annealing,  (b)  after  annealing  at  400°C  and  (c)  after  annealing  at  460°C.  K  =  -  1.6  V  for  all 
three.  The  area  is  about  150  x  150  for  (a)  and  (b),  and  about  43  x  43  for  (c),  respectively. 


272 


Xiang-Dong  Wang  ei  al.  /  Atomic  hydrogen  chemisorption  on  Si(100)(2  x  }) 


4b,  the  surface  is  characterized  by  a  high  density 
of  dark  defects  and  some  bright  dots  which  made 
the  surface  look  highly  disordered.  However,  the 
direction  of  the  dimer  chain  can  still  been  identi¬ 
fied,  indicating  that  it  is  of  the  (2x1)  reconstruc¬ 
tion.  We  attribute  this  change  to  desorption  of 
the  hydrogen  in  the  ^j-state,  while  the  monohy¬ 
dride  phase  still  exists.  There  exists  controversy 
about  the  way  of  Hj  desorption  from  the  dihy¬ 
dride  state.  Schulze  et  al.  [11]  and  Masson  et  al. 
[28]  had  suggested  that  the  desorbing  H2  was 
formed  by  the  association  of  the  two  hydrogen 
atoms  in  the  same  SiH2  species  in  order  to  ex¬ 
plain  the  suggested  first-order  desorption  kinet¬ 
ics.  In  contrast,  theoretical  work  [29]  assumed  the 
desorption  of  H2  from  the  dihydride  species  on 
Sid  00)  in  terms  of  the  recombinationary  desorp¬ 
tion  of  two  hydrogen  atoms  from  two  adjacent 
silicon  hydride  species.  In  this  case,  the  transition 
into  the  monohydride  phase  takes  place  sponta¬ 
neously  with  the  combination  of  hydrogen  mole¬ 
cules.  This  model  supported  the  second-order 
desorption  kinetics  obtained  from  the  isothermal 
desorption  experiment  [10].  However,  as  can  be 
seen  in  fig.  4b,  the  surface  is  highly  disordered 
due  to  the  interruption  of  the  dimer  chain  by 
dark  areas  which  are  attributed  to  the  monohy¬ 


dride  phase  from  our  previous  analysis.  This  sug¬ 
gests  that  H2  molecules  desorb  directly  from  the 
dihydride  phase  without  shifting  to  the  monohy¬ 
dride  phase.  Because  of  the  complicated  pro¬ 
cesses  involved  in  the  adsorption  and  desorption 
of  hydrogen,  though  it  might  be  difficult  to  rule 
out  unambiguously  one  of  the  desorption  models 
suggested  above,  we  suggest,  from  the  STM  ob¬ 
servation,  that  at  least  a  large  amount  of  hydro¬ 
gen  in  the  dihydride  phase  desorbs  directly  from 
the  surface  without  transiting  into  the  monohy¬ 
dride  phase.  Thus  our  STM  result  supports  the 
assumption  that  the  desorbing  H2  molecules  are 
formed  by  the  tw'o  hydrogen  atoms  in  the  same 
SiH2  species  instead  of  two  adjacent  SiHj  species. 
This  is  in  contradiction  with  TPD  data  [10,30], 
however,  TPD  data  did  not  give  an  unambiguous 
answer  to  the  kinetics  of  H2  desorption  in  the 
dihydride  phase  as  no  TPD  work  gave  the  abso¬ 
lute  coverage  change  before  and  after  desorption 
of  the  dihydride  phase.  In  fact,  the  desorption 
kinetics  of  H2  from  the  Si(lOO)  surface  has  not 
been  well  understood  yet,  the  previous  analysis  of 
TPD  data  was  merely  based  on  the  assumption 
that  the  dihydride  phase  transferred  to  monohy¬ 
dride  phase  by  desorbing  H2  from  two  dihydride 
species.  As  to  our  knowledge,  only  Oura  et  al. 


■<0i  i> 


(c) 


SuDlayer  Si  atom, 


OUnouckleCly 

aimerisea  atom, 


Hyflrogen  atom, 


•  O 


Upper  and  lower  atoms  in 
a  slightly  Puciclea  dimer. 


•  O 


Upper  and  lower  atoms  in 
a  seve'-ely  buckled  dimer. 


Fig.  5.  Schematic  of  the  adsorption  geometry  of  hydrogen  atoms,  (a)  initial  stage  of  adsorption  (proposed  in  this  work),  (h)  (2x1) 
monohydride  phase  (ref.  12]),  and  (c)  (1x1)  dihydride  phase  (ref.  [2]). 
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their  work  [31]  shows  the  absolute  coverage 
change  of  the  adsorbate  from  1.85  ML  of  hydro¬ 
gen  saturated  (1  x  1)  phase  to  1  ML  of  the  (2  x  1) 
monohydride  phase.  However,  this  result  dose 
not  have  enough  accuracy  to  convince  people 
since  an  ideally  saturated  (1x1)  phase  should  be 
2  ML  instead  of  1.85  ML.  Furthermore,  no  other 
work  has  ever  been  published  to  support  this 
conclusion. 

Upon  further  annealing  of  the  sample  at  460°C 
for  1  min,  the  surface  is  as  shown  in  fig.  4c,  which 
is  dominated  by  the  (2x1)  structure  with  or¬ 
dered  dimer  chains,  although  the  density  of  de¬ 
fects  is  still  higher  than  that  on  the  clean  surface. 
A  novel  feature  in  fig.  4c  is  the  formation  of  the 
adatom  lines  on  the  substrate  terraces,  while  the 
number  of  bright  dots  decreases  considerably, 
indicating  the  re-formation  of  the  clean  Si(100K2 
X  1)  surface.  By  correlating  the  disappearance  of 
the  bright  dots  and  the  formation  of  the  adatom 
lines,  we  attribute  these  adatom  lines  to  Si  atoms 
that  formed  by  desorption  of  hydrogen  from  SiHj 
and  small  amounts  of  other  species.  The  mobility 
of  dissociated  Si  atoms  is  high  enough  so  that  the 
Si  atoms  can  combine  with  each  other  with  high 
probability',  forming  epitaxial  dimer  lines  perpen¬ 
dicular  to  the  direction  of  the  dimer  lines  of  the 
substrate.  Further  annealing  of  the  sample  at 
500°C  results  in  no  obvious  changes  in  the  surface 
structure,  indicating  that  the  critical  temperature 
for  Si  overlayer  formation  is  between  400°C  and 
460°C. 


4.  Summary 

Our  STM  investigations  of  H  adsorption  on 
the  Sid  00X2  x  1)  surface  can  be  summarized  as 
follows  (see  fig.  5  for  illustration). 

(1)  At  low  coverages  (below  10%  of  a  mono- 
layer),  single  hydrogen  atoms  reside  on  the  top  of 
the  dimerized  surface  Si  atoms,  imaged  as  bright 
spots.  Those  H  atoms  also  induce  the  buckling  of 
the  dimers  at  both  sides  of  the  adsorbed  hydro¬ 
gen  atom. 

(2)  With  increasing  coverage,  the  density  of 
hydrogen  atoms  which  adsorb  in  pairs  increases. 
The  monohydride  phase  is  imaged  dark  com¬ 


pared  with  unreacted  Si  dimers  because  of  the 
reduction  of  the  density’  of  states  near  the  Fermi 
level. 

(3)  Our  data  suggest  that  the  reacted  Si  atoms, 
mainly  SiHj,  are  imaged  as  bright  dots,  spreading 
randomly  on  the  dimer  chain.  At  the  monohy¬ 
dride  phase,  the  surface  etching  is  slight,  while 
upon  the  formation  of  the  drhydride  phase,  the 
corrosion  of  the  surface  becomes  significant.  It 
seems  that  the  etching  originates  and  expands 
randomly  without  any  preferred  sites. 

(4)  Upon  annealing  a  surface  having  a  mixture 
of  the  (1x1)  and  (2x1)  phases  at  different 
temperatures,  the  changes  in  the  surface  struc¬ 
ture  were  observed.  Upon  annealing  at  400°C,  the 
dihydride  phase  disappeared  while  the  monohy¬ 
dride  phase  still  existed.  Our  results  support  the 
assumption  that  the  desorbing  H,  molecules  are 
formed  by  combination  of  two  hydrogen  atoms  in 
the  same  SiHj  species.  Upon  annealing  the  sam¬ 
ple  above  460°C,  the  (2x1)  monohydride  phase 
disappeared,  the  surface  is  dominated  by  the 
Si(2  X  1)  dimers  with  overlayer  stripes  formed  by 
the  recombination  of  Si  atoms  through  the  disso¬ 
ciation  of  SiHj  species.  Sponsor  [32]. 
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